SOLUTION OF QUARTIC EQUATIONS 

Program C2-35.0 
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REVISION NOTICE 

This publication replaces previous descriptions of "Solution of Quartic 
Equations," program C2-35.0. Descriptions have been clarified and pro- 
gram references changed to current nomenclature. 



FUNCTION 

"Solution of Quartic Equations" is used to determine all roots (real 
and imaginary) of the 4th Order Polynomial 

agz^ + a^z^ + a£z + a3Z + &4 = 



INPUT 

1- Store the program in Lo. 

2. Arrange aQ = 1. 

3. Store data in Lo + 0600. 

4. Store flags, counters, and constants as indicated on "DATA LOAD 
SHEET," page 4. 

5. Until the program is revised it will be necessary to store zeros 
in LO + 0605 to Lo + 0635, inclusive. 



OUTPUT 

Set tab stops for a 4 column printout. The printout will be in the 
form A = Bi . 
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SOLUTION OF QUARTIC EQUATIONS 



PROCEDURE 

Transfer to beginning of program and start with the MANUAL INPUT lever 
on the typewriter depressed. The light on the typewriter will then 
come on and the operator can feed the prepared data through the reader 
by lifting the typewriter MANUAL INPUT lever. The computer will store 
the data, make a carriage return, and then start to compute the roots. 

After the roots have been printed out it will be necessary to depress 
the START button on the computer and the following operations will be 
executed: 

1. The given coefficients will be printed out. 

2. Check coefficients will be computed from the roots. 

3. The check coefficients will be printed out. 

If the coefficients do not compare closely it is necessary to change 
the flag in Lo + 0640 to a higher number and repeat the entire process. 



EXAMPLE I 

In the event a program should stop due to an excessive exponent it is 
necessary to scale the coefficients of the equation. 

Use the following method: 

Given 

f(x) = x n + a^"- 1 + a 2 x n " 2 + a 3 x n - 3 + + a n 

Substitute 

10 m y = x 
Then 

n-1 n-2 n-3 
a iY + a 5 y a^y 
f (y) = y n + _i^ -K=- -^L_ + + _^_ 

10 m 10 2m 10 3m ^nm" 



2 of 4 



SOLUTION OF QUARTIC EQUATIONS 



EXAMPLE I (Cont. ) 



In terms of its roots f(y) can be expressed 

f( y ) = (y - r x ) (y - r 3 ) ( y _ r ) 



but. 



and 



x 



r f(x) = -T- C(X = "^X* " 10 m r 2 )(x - 10 m r )...(x - 10 m r n )] 



io ra 10' 



m 



From the above it can be easily seen that the roots of f(y) are l/io' 
of the roots of f(x). 



EXAMPLE II 

Cubic equations have been set up for this program to solve by multi- 
plying the cubic by (x - r), where (x - r) is a factor that is selected 
at random. 

Since f(x) - (x - r) f(x), it can be seen that by eliminating (x - r) 
from the roots of the resultant quartic we have the roots of the cubic. 



NOTE 



This entire subroutine is coded in floating point. The Floating Point 
Interpretive System 1 , program Hl-24.0, must be in memory and the first 
two instructions of this subroutine must be addressed to enter it. 
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SOLUTION OF QUARTIC EQUATIONS 
DATA LOAD SHEET 



NOTES 
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